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Abstract. We present the results of a systematic search for X- 
ray flares on young stars observed during ROSAT PSPC obser- 
vations of the Taurus-Auriga-Perseus sky region. All pointed 
PSPC observations currently available from the ROSAT Public 
Data Archive with known pre-main sequence T Tauri Stars or 
young Pleiads or Hyads in the field of view are analyzed. A 
study of the activity of late-type stars of different ages provides 
information on the evolution of their coronal activity, which 
may be linked to their angular momentum. 

We develop a criterion for the detection of flares based 
on the shape of the X-ray lightcurve. Applying our detection 
method to all 104 PSPC pointings from the archive we find 
52 flares. Among them 15 are detected on T Tauri Stars, 20 
on Pleiads, and 17 on Hyads. Only the 38 events which can 
definitely be attributed to late-type stars (i.e. stars of spectral 
type G and later) are considered in the statistical analysis of the 
properties of flaring stars. We investigate the influence of stellar 
parameters such as age, rotation and multiplicity on individual 
flare parameters and flare frequency. 

From the total exposure time falling to the share of each 
sample and the duration of the individual flares we compute a 
flare rate. We take into account that the detection sensitivity for 
large X-ray flares depends on the S/N and hence on the stellar 
distance. The values we derive for the flare rates are 0.86 ± 
0.16% for T Tauri Stars, 0.67 ± 0.13% for Pleiads and 0.32 ± 
0.17% for Hyads. The flare rate of classical T Tauri Stars may 
be somewhat higher than that of weak-line T Tauri Stars (Fc — 
1.09 ± 0.39% versus = 0.65 ± 0.16%). 

Hardness ratios are used to track the heating that takes place 
during stellar flares. Hardness ratios are evaluated for three dis- 
tinct phases of the flare: the rise, the decay, and the quiescent 
(pre- and post-flare) stage. In most cases the hardness increases 
during the flares as compared to the quiescent state. During 
both quiescence and flare phase TTSs display the largest hard- 
ness ratios, and the Hyades stars show the softest spectrum. 

Key words: stars: flare - X-rays: stars - stars: late-type 



1. Introduction 

The Taurus-Auriga-Perseus region offers the opportunity to 
study the X-ray emission of young stars at several evolution- 
ary stages. The youngest stars observed by ROSAT in this por- 
tion of the sky are the T Tauri Stars (TTSs) of the Taurus- 
Auriga and Perseus star forming regions, late-type pre-main 
sequence (PMS) stars of M < 3Mq with an estimated age 
of 10^ — 10^ yrs. Two young star clusters, the Pleiades and 
Hyades, are also located in this region of the sky at age of 
10^ yrs and 6 10* yrs, respectively. They consist mostly of 
zero-age main-sequence (ZAMS) stars, except for some higher 
mass post-main sequence stars and brown dwarfs, which are 
not studied here. 

From the early observations by the Einstein satellite it was 
concluded that the X-ray emission of young stars arises in an 
optically thin, hot plasma at temperatures above 10^ K (Feigel- 
son & DeCampli 1981). The emission region has been associ- 
ated with the stellar corona where the X-rays are produced — 
more or less analogous to the solar X-ray emission — through 
a stellar a-fi-dynamo. The dynamo is driven by the combina- 
tion of rotation and convective motions. Correlations between 
the X-ray emission of late-type stars and the stellar rotation 
support the notion that dynamo-generated magnetic fields are 
responsible for heating the coronae ( ^allavicini et al. 1981 ). 
But successful direct measurements of the magnetic fields of 
TTSs have been performed only recently (see e.g. Guenther 
et al. 1999). The details of the heating mechanism are still not 
well understood. 

The correlation between stellar rotation and X-ray emis- 
sion of late-type stars suggests that the rotational evolution of 
young stars determines the development of stellar activity. The 
rotational evolution of low-mass PMS stars partly depends on 
the circumstellar environment. While classical TTSs (hereafter 
cTTSs) are surrounded by a circumstellar disk, infeiTed from 



IR dust emission (Bertout et al. 1988, Strom et al. 1989, and 



Beckwith et al. 1990) and more recently from direct imaging 



(e.g. McCaughrean & O'Defl 1996), weak-line TTSs (wTTSs) 
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lack such a disk, or at least the disk is not optically thick. Ow- 
ing to contraction wTTSs spin up as they approach the main 
sequence. For cTTSs, on the other ha nd, coupling betwee n the 
disk and the star may prevent spin-up ( Bouvieretal. 1993 ). The 



2 



B. Stelzer, R. Neuhauser & V. Hambaryan: X-ray flares in Taurus-Auriga-Perseus 



period observed on the ZAMS depends on the time the star has 
spent in the cTTS phase. After the main-sequence is reached, 



the rotation rate decreases again (see Bouvier et al. 1997c ). As 
a consequence of their slower rotation, stars on the ZAMS and 
main sequence (MS) should on average show less X-ray activ- 
ity than PMS stai's. 

Earlier investigations of X-ray observations of young late- 
type stars were mostly concerned with the quiescent emission 
(see l^euhauser et al. 1995 , Stauffer et al. 1994 , Gagne et al. 



1995, Hodgkinetal. 1995 



Micelaetal. 1996 1999, Pye et al 



1994, and Stern et al. 1994). In contrast to these studies we fo- 
cus on the occurrence of X-ray flares. Furthermore we discuss 
a larger sample than most of the previous studies by using all 
currently available observations from the ROSAT Public Data 
Archive that contain any TTS, Pleiad or Hyad in the field of 
view. 

X-ray flares may be used as a diagnostics of stellar activ- 
ity. They are thought to originate in magnetic loops. In con- 
trast to findings from quasi-static loop modeling, the only di- 
rect determination of the size of a flaring region (Schmitt & 
Favata 1999) shows that the emitting region is very compact. 
In the loops which confine the coronal plasma magnetic recon- 
nection suddenly frees large amounts of energy which is dis- 
sipated into heat and thus leads to a temporary enhancement 
of the X-ray emission. The decay of the lightcurve is accom- 
panied by a corresponding (exponential) decay of the temper- 
ature and emission measure, which are obtained from one- or 
two-temperature spectral models for an optically thin, thermal 
plasma ( [Raymond & Smith 1977] , ^ewe et al. 1985] , 1986). 



The most powerful X-ray flares have been observed on the 
youngest objects, notably a flare on the infrared Class I pro- 
tostar YLW 15 in p Oph which has been presented by Grosso 
et al. (1997). X-ray flares on TTSs observed so far (see Mont- 



merle et al. 1983, 


Preibisch et al. 1993 




Strom & Strom 1994, 


Preibisch et al. 1995 


, Gagnc et al. 1995 




Skinner et al. 1997, 


Tsuboi et al. 1998 


) exceed the maximum emission observed 



from solar flares by a factor of 10 and more. Some extreme 
events have shown X-ray luminosities of ~ Lx = 10^"^ erg/s. 
Although some of the strongest X-ray flares ever observed were 
detected on TTSs to date no systematic search for TTS flares 
was undertaken. 

This paper is devoted to a study of the relation between 
X-ray flare activity and other stellar parameters, such as age, 
rotation rate, and multiplicity. For this purpose we perform 
a statistical investigation of ROSAT observations. We develop 
a method for the flare detection based on our conception of 
the typical shape of a flare lightcurve, where the term 'typi- 
cal shape' refers to the characteristics of the X-ray lightcurve 
described above, i.e. a significant rise and subsequent decay 
of the lightcurve to the previous emission level. The database 
and source detection is described in Sect. ^. In Sect. |^ we de- 
scribe how the lightcurves are generated. Our flare detection 
algorithm is explained in Sect. where we also present all 
flare parameters derived from the X-ray lightcurves. Then we 
describe the influence of observational restrictions on the data 
analysis and how the related biases can be overcome (Sect. 1). 



In Sect. ^ we compare the flare characteristics of different sam- 
ples of flaring stars selected by their age, rotation rate, and mul- 
tiplicity. We present luminosity functions for TTSs, Pleiads, 
and Hyads during flare and quiescence. Luminosity functions 
of the non-active state of these stars have been presented before 



(see e.g. Pyeetal. 1994, Hodgkinetal. 1995, Neuhauser 1997) 



and some of the flares discussed here have been discussed by 
pagnc et al. (1995)| , |s"trom & Strom (1994)| , and Preibisch et al. 
(1993). However, this is the first statistical evaluation of flare 
luminosities. Flare rates comparing stellar subgroups with dif- 
ferent properties (such as age, v sin i, and stellar multiplicity) 
are compiled in Sect. ^ Because of lack of sufficient statistics 
for a detailed spectral analysis, hardness ratios are used to de- 
scribe the spectral properties of the flares. In Sect. || we present 
the observed relations between hardness ratios measured dur- 
ing different activity phases and between hardness and X-ray 
luminosity. Finally, we discuss and summarize our results in 
Sect. I and Sect. 

2. Database and data reduction 

In this section we introduce the stellar sample and explain the 
analysis of the raw data. Details about our member ship lists 
for TTSs, Pleiads, and Hyads are given below (Sect. 2A). We 
have retrieved all pointed ROSAT PSPC observations from the 
archive that contain at least one of the stars from these lists in 
their field. The observations are listed in Table ^ After per- 
forming source detection on all of these pointings, we have 
cross correlated the membership lists with the detected X-ray 
sources and identified individual TTSs, Pleiads, and Hyads in 
the X-ray image. The process of source detection and identifi- 
cation is described in Sect. 2.2 



2.1. The stellar sample 

The analysis presented here is confined to the Taurus-Auriga- 
Perseus region. This portion of the sky includes the Taurus- 
Auriga complex, the MBM 12 cloud, and the Perseus molecu- 
lar clouds with the reflection nebula NGC 1333 and the young 
cluster 1C348. Two open clusters containing mostly ZAMS 
stars are found nearby the above mentioned star forming re- 
gions, the Pleiades and the Hyades. The choice of this specific 
sky region thus enables us to compare the X-ray emission of 
young stars at different ages. 

Our sample of low-mass PMS stars in and around Taurus 
consists of all TTSs which are either on or very close to the Tau- 
rus star forming clouds or off the clouds at locations where they 
can still be linked with the Taurus clouds (see e.g. Neuhauser 
et al. 1997). We restrict our Taurus sample to objects between 
0^2000 = 2/i and 5h and 62000 = —10° and 40°. The TTSs in 
Taurus comprise those listed in the Herbig-Bell catalog (Herbig 



& Bell 1988; HBC), in Neuhauser et al. (1995) or in Kenyon 



& Hartmann (1995). In addition we include TTSs newly iden- 
tified eit her as counterparts to p reviously unidentified R OSAT 
sources ( [Strom & Strom 1994|, [Wichmann et al. 1996t Mag- 
azzii et al. 1997, [Neuhauser et al. 1997| , [Zickgraf et al. 1998 , 
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Fig. 1. Sky map of the Tau-Aur-Per region showing the point- 
ing positions of the ROSAT PSPC observations. On display 
are all PSPC pointings selected from the ROSAT Public Data 
Archive which contain at least one TTS, Pleiad, or Hyad from 
the membership lists given in the text. In addition all stars from 
the membership lists are plotted. 



Li & Hu 1998|, [Briceno et al. 1998 ) or by other means (T orres 



etal. 1995, pppenheimeretal. 1994 



Briceno etal. 1998, Reid 



& Hawley 1999, Gizis & Reid 1999). TTSs in the molecular 



cloud MBM 12 (see [Hearty et al. 2000[ ) are also in the exam- 
ined sky region. 

In addition to TTSs from Tau-Aur we include those from 
the Perseus molecular cloud complex, mainly IC348 and 
NGC 1333, in our analysis. Our list of TTSs in IC348 com- 
prises X-ray detections identified with Ha emission stars or 
with proper motion members (Tables 4 and 5 in Preibisch et al. 
1996), and emission line stars from Herbig (1998), and Luh- 
man (1999). TTS members of NGC 1333 ai'e listed in Preibisch 
(1997). 

All objects with low lithium strength are excluded, because 
it is dubious whether they are young. We accept only those 
objects as PMS stars which show more lithium than Pleiades 
stars of the same spectral type, i.e. we exclude all those with 
WA(Li) lower than 0.2 A for F- and G-type stars and lower than 
0.3 A for K-type stars. When applying this criterion, we always 
use the spectrum with the best resolution and best S/N, i.e. the 
high-resolution spectra from Wichmann et al. (in preparation). 
If no high-resolution sp ectra are available, we use th e medium- 
resolution spect ra from ^artm & Magazzu (1999) , Neuhauser 
et al. (1997), or |Magazzii et al. (1997)| . 

Members of the Pleiades and Hyades clusters are se- 
lected from the Open Cluster database compiled by C. 



Prosser and collegues (available at http://cfa-www.harvard.edu/ 



^stauffer/opencl/index.html). The tables collected in Prosser's 
database provide a summary of membership classification 
based on different methods, such as photometry, spectra, ra- 
dial velocity and Ha emission. In addition a final membership 
determination is given which we use to define our membership 
Usts. 



Fig. |l] shows a sky map with the positions of the selected 
PSPC observations. 



2.2. Source detection and identification 

Source detection is performed on all observations given in Ta- 
ble |l] and shown in Fig. [l| using a combined local and map 
source d etection algorithm b ased on a maximum likelihood 
method ( [Cruddace et al. 1988| ). All detections with ML> 7.4 
(corresponding to ^ 3.5 Gaussian a determined as best choice 



by S^euhauser et al. 1995) are written to a source list, which is 
subsequently cross-correlated with the membership lists intro- 
duced above. The maximum distance A between optical and X- 
ray position to be allowed in this identification process depends 
on the off-axis angle of the source because the positional accu- 
racy of the PSPC is worse at larger distances from the center 
due to broader point spread function (PSF). From distributions 
of the normalized cumulative number of identifications ver- 
sus offset A for different off-axis ranges we have determined 
the optim um cross-correlation ra dius for all detector positions 
(similar to 'Neuhauser et al. 1995 ). A detailed description of this 
process together with a table providing A for different off-axis 
ranges will be given in Stelzer et al. (in preparation). 

Observations which are characterized by strong back- 
ground variations (200020, 200008-0, and 200442), as well as 
observations consisting of two very short intervals separated by 
a gap of > 100 h (200068-0, 200914) ai'e omitted from the flare 
detection and flare analysis. Furthermore, we neglect observa- 
tions with a total duration of less than 1000 s. All observations 
which have been ignored in the analysis presented in this paper 
are marked with an asterisk in Table |l| 



3. Lightcurves 

Using the arrival time information of the photons counted 
within a pre-defined source circle lightcurves are generated for 
each of the ROSAT sources that have been identified with a 
TTS, Pleiad or Hyad from the membership lists. 

For the source extraction radius we have used the 99% 
quantile of the Point Spread Function (PSF) at 1 keV, i.e. the 
radius containing 99% of the 1 keV photons at the respective 
off-axis angle. In contrast to the standard EXSAS source ra- 
dius of 2.5 FWHM, which becomes unreasonably large for off- 
axis sources due to the extended wings of the PSF, this choice 
of extraction radius limits the source size. Close to the detector 
center, some bright sources slightly overshine the nominal 99% 
quantile of the PSF probably due to small deviations from the 
assumed 1 keV spectrum. We have therefore checked all im- 
ages for such bright sources and determined a larger source ra- 
dius for these cases based on visual inspection. In crowded re- 
gions, where the PSF of several sources overlap, the measured 
counts are upper limits to the actual emission of the sources. 
None of the overlapping sources showed a flare, however, such 
that no further attention is drawn to the overestimation of the 
count rate in these cases. 
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Table 1. Complete list of ROSAT PSPC pointings available 
from the ROSAT V\x\i\\c Data Archive in October 1998, which 
include at least one TTS, Pleiad or Hyad from our member- 
ship lists (see text) in the field of view. All observations not 
marked by an asterisk have been analyzed in this work. The to- 
tal exposure time is given in Column 3. The remaining columns 
give the number of detected and undetected TTSs, Pleiads and 
Hyads (D - detection, N - non-detection). If more than one 
star was identified with an X-ray source all identifications are 
listed. For multiple systems each component is counted as one 
entry. 



No. 


ROSAT 


exp (s) 


TTS 


Pleids 


Hyads 




Pointing ID 




D/ N 


D/ 


N 


D/ 


N 


1 


180185p 


8897 


16/ 7 


0/ 





0/ 





2 


200001p-l 


1520 


0/ 2 


0/ 





0/ 





3 


200001p-0 


4486 


23/12 


0/ 





0/ 





4 


200001p-l 


25591 


26/ 9 


0/ 





0/ 





5 


200008p-2* 


696 


0/ 6 


2/162 


0/ 





6 


200008p^* 


121 


0/ 5 


/ 220 


0/ 





7 


200008p-5* 


722 


0/ 4 


1/110 


0/ 


2 


8 


200008p-0* 


5936 


21 4 


86 / 165 


1/ 





9 


200008p-2 
200020p* 


7049 


4/ 3 


88/161 


1/ 





10 


39879 


0/ 2 


0/ 





13/ 


8 


11 


200068p-l 

200068p-0* 


1307 


3/ 2 


49 / 184 


0/ 





12 


12849 


3/ 2 


93/ 


155 


1 / 





13 


200068p-l 


27071 


3/ 2 


101 / 145 


1 / 





14 


200073p 


2376 


0/ 1 


0/ 





6/ 


3 


15 


200082p-r 


814 


0/ 


0/ 





0/ 


2 


16 


200083p 


2799 


0/ 


0/ 





16/ 


5 


17 


200107p-l 


27692 


0/ 


0/ 





21 


1 


18 


200107p-0 


3923 


0/ 


0/ 





21 


1 


19 


200402p 


10469 


3/ 


0/ 





0/ 





20 


200441p 


10987 


0/ 3 


0/ 





15/ 


8 


21 


200442p 


19948 


21 


0/ 





6/ 


4 


22 


200443p 


20074 


15/ 8 


0/ 





12/ 


2 


23 


200444p 


14593 


0/ 1 


0/ 





21 


5 


24 


200547p 


28359 


0/ 3 


0/ 





0/ 





25 


200553p 


10961 


0/ 


0/ 





8/ 


3 


26 


200556p 


22456 


5/ 


67/ 


81 


3/ 





27 


200557p 


27648 


7/ 


94/ 


96 


0/ 





28 


200576p 
200677p* 


1521 


0/ 


0/ 





12/ 


2 


29 


650 


3/13 


0/ 





3/ 





30 


200694p 


1987 


17/13 


0/ 





1 / 


1 


31 


200694p-l 


5395 


16/ 12 


0/ 





21 





32 


200775p 


4096 


0/ 2 


0/ 





3/ 


9 


33 


200776p 


22995 


0/ 1 


0/ 





10/ 


5 


34 


200777p 


16296 


0/ 1 


0/ 





13/ 


8 


35 


200778p 


1900 


0/ 1 


0/ 





15/ 


11 


36 


20091 Ip 


17460 


0/ 


0/ 





7/ 


4 


37 


200911p-l 


13755 


0/ 


0/ 





7/ 


4 


38 


200912p 


1656 


0/ 2 


0/ 





3/ 


6 


39 


200912p-l 


23710 


1/ 1 


0/ 





5/ 


4 


40 


200913p 
200914p* 


25341 


21 


0/ 





7/ 


2 


41 


4098 


0/ 


0/ 





1 / 


1 


42 


200915p 


3504 


0/ 


0/ 





21 


2 


43 


200942p 


7367 


0/ 


0/ 





M 


6 


44 


200945p 


4099 


0/ 


0/ 





11 


4 


45 


200949p 


6098 


17/ 3 


0/ 





0/ 





46 


200980p 


10649 


0/ 


0/ 





3/ 


3 


47 


200980p-l 


4014 


0/ 


0/ 





3/ 


3 


48 


200981p 


4491 


21 1 


0/ 





8/ 


7 


49 


200982p 


7724 


0/ 


0/ 





10/ 


3 


50 


201012p 


7492 


21 1 


0/ 





0/ 





51 


20101 3p 


4826 


3/ 2 


0/ 





0/ 





52 


201013p-l 


4484 


21 4 


0/ 





0/ 





53 


201014p 


9910 


11 1 


0/ 





1/ 


1 


54 


201015p 


10058 


0/ 1 


0/ 





1/ 


1 


55 


201016p 


10576 


6/ 6 


0/ 





0/ 





56 


201017p 


8452 


7/ 9 


0/ 





0/ 





57 


201023p 


3058 


0/ 


0/ 





0/ 


3 


58 


201025p 


5448 


11/17 


0/ 





0/ 





59 


201097p 


10287 


0/ 


0/ 





1 / 





60 


201278p 


1358 


4/ 


0/ 





0/ 





61 


201278p-l 


4028 


4/ 


0/ 





0/ 





62 


201305p 


23643 


94/96 


0/ 





0/ 





63 


201312p 


2800 


8/26 


0/ 





0/ 





64 


201313p 


4027 


22/12 


0/ 





11 


6 


65 


2013 14p 


2691 


0/ 


0/ 





Yll 


4 


66 


201314p-l 


1397 


0/ 


0/ 





91 


6 



Table 1. continued 



No. 


KUoAl 


exp (s) 


T 




p 


H 




Pointing ID 




D/N 


D/N 


D/N 


67 


201315p'' 


645 


4/ 


1 


0/ 


1/9 


68 


201315p-l 


1544 


5/ 


0/ 


4/7 


69 


201315p-2 


2298 


5/ 


0/ 


4/ 6 


70 


20l3l6p 


4156 


1/ 


1 


0/0 


3/3 


71 


201317p 


1731 


0/ 


0/ 


1/ 1 


72 


201319p 


1811 


4/ 


0/ 


3/ 1 


73 


201368p 


16594 


0/ 


1 


0/ 


17/9 


74 


201369p 


15538 


0/ 


0/ 


24/ 4 


75 


201370p 


4975 


0/ 


1 


0/0 


5/4 


76 


201370p-l 


13652 


1/ 


0/0 


5/2 


11 


201484p 


7673 


1/10 


0/0 


0/0 


78 


201485p 


2429 


0/ 


3 


0/ 


0/ 


79 


201485p-l 


1904 


1/ 2 


0/0 


0/0 


80 


201504p 


109534 


2/ 


0/0 


4/0 


81 


201519p 


5889 


2/ 


0/ 


0/ 


82 


201532p 


10183 


9/ 6 


0/0 


1/0 


83 


201533p 


10799 


4/ 6 


0/0 


0/0 


84 


201534p 


6213 


0/ 


1 


0/ 


3/ 1 


85 


2U15ysp 


5652 


7/14 


0/ 


0/ 


86 


201599p 


6175 


6/15 


0/0 


0/0 


87 


201600p 


5747 


6/15 


0/0 


0/0 


88 


201601p 


5817 


5/16 


0/0 


0/0 


89 


201602p 


5582 


8/13 


0/ 


0/ 


90 


201747p 


19600 


0/ 


1 


0/0 


16/2 


91 


201748p 


16999 


0/ 





0/0 


2/3 


92 


2{)1749p* 


924 


0/ 





0/ 


3/ 1 


93 


201749p-l 


1540 


0/ 





0/ 


3/ 1 


94 


300178p 


710 


0/ 


7 


0/0 


0/0 


95 


4003 12p 


10735 


0/ 





0/ 


1/0 


96 


700044p 


461 1 


5/ 





0/ 


61 4 


97 


700063p 


1701 


0/ 





0/0 


1/0 


98 


700825p 


1435 


0/ 





0/0 


0/2 


99 


700825p-l 


15744 


0/ 





0/ 


1/ 1 


100 


7009 13p 


2096 


0/ 





0/0 


1/ 


101 


7009 1 6p 


7161 


0/ 





0/0 


1/ 


102 


7009 1 9p 


2004 


0/ 





0/0 


1/0 


103 


/uuy^jp 


25 16 


0/ 





0/ 


1/ u 


104 


701055p 


9039 


0/ 





0/0 


0/0 


105 


701253p 


5417 


0/ 





0/0 


0/2 


106 


80005 lp-0 


1470 


1/ 





0/0 


0/0 


107 


800051p-l 


3347 


1/ 





0/0 


0/0 


108 


80005 lp-2 


3303 


1/ 





0/0 


0/0 


109 


800083p 


10220 


1/ 





0/0 


0/0 


110 


800104p 


7524 


1/ 





0/0 


0/0 


111 


800193p 


21969 


0/ 





0/0 


3/3 


112 


9001 3 8p 


24899 


6/ 





0/0 


0/0 


113 


900154p 


28791 


0/ 


2 


0/0 


0/0 


114 


900193p 


9128 


3/22 


0/0 


0/0 


115 


900353p 


7718 


18/11 


0/0 


12/4 


116 


900371p 


4227 


1/ 


0/0 


2/0 


117 


900371p-l 


7718 


0/ 


0/0 


2/0 



The events measured within the circular source region are 
binned into 400 s intervals. Since the typical duration of a flare 
is less than one hour, significantly longer integration times 
would lead to a loss of information about the structure of the 
Ughtcurve, while for shorter bin lengths the Ughtcurves are 
dominated by the low statistics. Furthermore, the choice of 
400 s integration time guarantees that no additional variability 
is introduced by the telescope motion (wobble). 

Due to the earth eclipses the data stream is interrupted at 
periodical time intervals. Depending on the phase used for the 
time integration, at the beginning and/or end of each data seg- 
ment the 400 s intervals are only partly exposed. For the flare 
detection only bins with full 400 s of exposure are used. To gain 
independence of the biiming we generate Ughtcurves with dif- 
ferent phasing of the 400 s intervals: First, in order to divide the 
given observing time into as many 400 s exposures as possible. 
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a lightcurve is binned in such a way that a new 400 s inter- 
val starts after each observation gap. Thus, data are lost only 
at the end of each data segment, because the last bin remains 
uncomplete. Secondly, lightcurves are built by simply splitting 
the total observing time into 400 s intervals beginning from the 
start of the observation regardless of data gaps. In this case data 
are rejected at the beginning and end of each data segment. 

The number of background counts falling in the source cir- 
cle is determined from the smoothed background image which 
is created by cutting out the detected sources and then perform- 
ing a spline fit to the resulting image. This method of back- 
ground acquisition is of advantage in crowded fields where 
an annulus around the source position - the most widely used 
method for estimating the background - likely is contaminated 
by other sources. The background count rate is found by di- 
viding the number of background counts in the source circle 
through the exposure time extracted from the standard ROSAT 
exposure map. To take account of possible time variations in 
the background count rate, the background is determined sepa- 
rately for each data segment and subtracted from the measured 
count rate in the respective data interval. (When referring to 
'data segments' we mean parts of the lightcurve that are sepa- 
rated from each other by gaps of at least 0.5 h.) 

4. Flare detection 

4.1. The method 

One of the major elements of a flare by customary definition is 
a significant increase in count rate, after which the initial level 
of intensity is reached again. Therefore, our flare detection is 
based on the deviation of the count rate from the (previously 
determined) mean quiescent level of the source. To ensure that 
the quiescent count rate contains no contribution from flares, 
in the first step, we determine mean count rates for all data 
segments of each lightcurve and define the quiescent level as 
the lowest mean measured in any of these data segments. 

We define a flare as an event which is characterized by two 
or more consecutive time bins that constitute a sequence of ei- 
ther rising or falling count rates, corresponding to rise and de- 
cay phase of the flare. In addition, to ensure the significance of 
our flare detections, we define the upper standard deviation of 
the quiescent level as a point of reference and require that (a) all 
bins which are part of the flare are characterized by count rates 
higher than this level, and that (b) the sum of the deviations of 
all these bins is more than 5 a from this level. A rise imme- 
diately followed by a decay is counted as one flare. Since the 
shape of a lightcurve is influenced to some degree by the bin- 
ning used, we accept only flares that are detected in lightcurves 
with both bin phasings (see Sect. 

Detections of more than one flare in a single lightcurve are 
possible. To estimate the contribution of each event properly, 
after detection the decay of the first flare in each lightcurve is 
modeled by an exponential function, and a new lightcurve is 
generated by subtracting the fit function from the data. Having 
removed the first flare, we search for further flares in the re- 
duced, 'flare-subtracted' lightcurve using the same criteria as 



before. This procedure is repeated until no additional flares are 
detected. 

Since many of the investigated sources are highly variable 
X-ray emitters on timescales shorter than resolvable by our 
method, the mean count rate used until now in some cases is 
not a good estimate for the quiescent emission. With the knowl- 
edge obtained about the times at which flares have occurred we 
therefore redetermine the quiescent count rate taking the mean 
from the remaining data after removal of all flare contributions. 
Using this new mean count rate we repeat the flare detection 
procedure. 

4.2. Flare Parameters 

With the detection procedure described in the previous sub- 
section we have found 52 flares. We have always identified 
the nearest optical position with the X-ray source. In one 
flare, however, two possible optical counterparts, DD Tau and 
CZTau, are closeby (at 6" and 24" respectively), so that we 
can not be sure which star flared. Fifteen events were ob- 
served on TTSs, 20 on Pleiads, and 17 on Hyads. On two TTSs 
(RXJ0437.5H-1851 and TTau) and two Hyads (VA334 and 
VB 141) two flares occurred in the same observation. VB 141 
showed a third event during a different /?(9&4r exposure. 

Hyades stars above 2 Mq, that have already evolved off the 
main-sequence, are not considered in the statistical analysis if 
they showed a flare. Brown dwarfs in the Hyades and Pleiades 
are not on the main-sequence per definition, but they are also 
too faint for X-ray detection (Neuhauser et al. 1999). Thus we 



discuss only the ZAMS from the Pleiades and Hyades. 

A complete list of all TTSs, Pleiads, and Hyads on which at 
least one flare was detected is given in Table ^. Column 1 gives 
the designation of the flaring star Column 2 is the distance es- 
timate used for the count-to-energy-conversion. For TTSs in 
Taurus-Auriga we adopt a value of 140 pc ( |31ias 1978 , Wich- 
mann et al. 1998), while the TTSs in MBM 12 are located at 



65 pc (Hearty et al. 2000), and those in Perseus are located at 
350pc (NGC 1333; ^erbig & Bell 1988[ ) and 300 pc (IC348; 
[Cernicharo et al. 1985 ). Pleiads are assumed to be at a distance 
of 116 pc, the value derived by Mermilliod et al. (1997). Fi- 



nally, we use the individual Hipparcos parallaxes for Hyades 
stars if available, and otherwise the mean value of 46 pc (Per- 
ryman et al. 1998). We give spectral type, v sini, multiplicity, 
and binary separation of the stars and their respective refer- 
ences in columns 3-9. For TTSs additional columns specify 
whether the star is a cTTS or a wTTS. 

The lightcurves of all flares have been analyzed by fitting 
an exponential to the decay of each flare. In Figs. ^, ||, and ^ 
the fit function and measured mean quiescent count rate are 
displayed together with the data. In Tables ^, ^ and || we give 
the result of the modeling. Column 1 and column 2 contain the 
stellar identification of the X-ray source and the ROSAT obser- 
vation request number (ROR). The mean quiescent count rate 
is given in column 3, the maximum count rate inferred from 
the exponential fit to the lightcurve in column 4, and the decay 
timescale Tdoc from the fit in column 5. For flares with poor 



6 



B. Stelzer, R. Neuhauser & V. Hambaryan: X-ray flares in Taurus-Auriga-Perseus 



t =^ is JD 2450507.999319 LkHa270(l«0185) 



is JD 2448321.338625 BPTau (200001-0) 




"Time [h| 

= is JD 24484S9.4804I5 HD283572 (200001-1) 



V 



Time |h| 

t ^ IS JD 2449046.083444 RXJ0437.5-I-1 8516(200913) 




Ois JD 2448321.338625 V410x-ray7 (200001-0) 



= is JD 2448489.480415 DDTau/CZTau (200001-1 ) 



0.30 

I 



5 ID 15 20 

1 = is JD 2449051.203142 LkCal9 (201278-1 1 




t - is JD 2448323. 1 24435 LI 55 1 -5 1(200443) 



t - is JD 2448859.926! 14 LHa92 (201305) 



= is JD 2448314.263116 RXJ0422. 1+1934 (700044) 




t^OisJD 2448314.263116 TTaii (700044) 



t = is JD 2448468.890683 RXJ025 5. 4+2005 (900138) 



= OisJD 2448679.881614 LkHa325 (900193) 




Fig. 2. PSPC lightcurves for flares on TTSs in Taurus-Auriga and Perseus detected by the procedure described in Sect. 4.1 



Identification of the X-ray source and ROSAT observation request number (in brackets) are given for each source. Binsize is 
400 s, 1 (T uncertainties. The dashed line represents the quiescent count rate. Solid lines are exponential fits to the data points 
belonging to the flare down to the quiescent emission. Background count rates are shown as upper limits when the background 
subtracted count rate is below zero. 



B. Stelzer, R. Neuhauser & V. Hambaryan: X-ray flares in Taurus-Auriga-Perseus 



7 



t = is JD 2448S6I.647213 hii 1 384(200008-2) 



= is JD 2448861.647213 hii2l47a+b(200008-2) 
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= is JD 2448491 .4066 1 9 hii303 (200068- 1 ) 
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Time |h| 

= is JD 244849 1 .4066 1 9 hcg307 (200068-1 ) 



Time I ill 

t - is JD 2448506.508770 hcg422 (200556) 



t = is JD 2448506.508770 hiil516 (200556) 



= Ois JD 2448505.510171 hcgl43 (200557) 



I = is JD 2448505.510171 sk702 (200557) 



iii-*-' it \ 



Ois JD 2448505.510171 hiil74 (2005571 






t 1 


(1 5 10 15 2(1 
Time |l>| 

t^OisJD 2448505.510171 hiil91 (200557) 


■ 1 


\ , 



I- 



OisJD 2448505.510171 hcg97 (200557) 




Fig. 3. PSPC lightcurves for flares on members of the Pleiades cluster detected by the procedure described in Sect 4. 1 . Identifi- 
cation of the X-ray source and /?(9&4r observation request number (in brackets) are given for each source. Binsize is 400 s, 1 a 
uncertainties. The dashed line represents the quiescent count rate. Solid lines are exponential fits to the data points belonging 
to the flare down to the quiescent emission. Background count rates are shown as upper limits when the background subtracted 
count rate is below zero. 



8 



B. Stelzer, R. Neuhauser & V. Hambaryan: X-ray flares in Taurus-Auriga-Perseus 



Ois JD 2448505.510171 hiil653 1200557) 




- is JD 2448505.5I0I7I Iiii345 (200557) 



OisJD 2448505.510171 hcglSl (200557) 



Time |h| 

t = is JD 2448505.510171 hii253 (200557) 




Tims I ill 

OisJD 2448505.510171 hii212 (200557) 



..-Hi i-ri---n--*i#i 



1 - 



]5 2U 



Fig. 3. continued 



data sampling we did not determine the errors of Tdcc- Col- 
umn 6 is the estimated rise time of the flare. Due to data gaps in 
most cases no reasonable estimate can be given. Luminosities 
are listed in columns 7 and 8: quiescent luminosity Lqui, and 
maximum luminosity during the flare L,nax- We assume that all 
stars in the system contribute the same level of X-ray emission 
during quiescence, but that only one component flares at any 
one time. Therefore, for all multiple stars the observed quies- 
cent count rate from column 3 has been divided by the number 
of components before the conversion to luminosity and energy. 

For the conversion from count rates to luminosities we 
have used th e mean ROSAT PSPC en ergy-conversion-factor 
(ECF) from [Neuhauser et al. (1995)| , i.e. ECF = 1.1 • 
10^^ cts cm^/crg and the distances given in Table ^ In or- 
der to eliminate uncertainties in the distance estimate we have 
computed ratios of luminosity (given in column 9). Here, ip 
denotes the luminosity emitted during the flare, i.e. Lp = 
Lmax — iqui- The total emitted energy during quiescence Eq^i 
(column 10) and during the flare alone E-p (column 11) are 
inferred from the integration of the lightcurve between imax 
and iniax + t. The last column gives the reference for flares 
that have been published previously. In the last two rows of Ta- 
bles || and I we have listed the mean and median for each of 
the given parameters, except Tiis which is not well constrained. 
The means and medians have been computed with the ASURV 
Kaplan-Meier estimator (see [Feigelson & Nelson 1985 ), tak- 
ing account of upper/lower limits. Lower limits of L,nax oc- 
cur when there is doubt about whether the maximum emission 
of the flare has been observed (due to a data gap near the ob- 
served maximum). Upper limits for Tdcc occur when the de- 
cay is not observed because of a data gap between maximum 



and post-flare quiescent count rate. Both luminosity and decay 
timescale determine the flare energy, but the limits of these two 
parameters carry opposite signs. We therefore consider all val- 
ues of E-p uncertain where Tdcc or imax are a limit (indicated 
by colons in Tables ^ Q and ^) and have not included them in 
the computation of the mean and median. 

5. Observational selection effects 

It is the purpose of this paper to compare the flare activity of 
different stars, and thus some attention has to be drawn to ob- 
servational selection effects. In this section, we will discuss 
how observational restrictions influence the search for flares. 
At several points during the data analysis, we are confronted 
with the problem of finding a representation of the data which 
is free from these biases. 

The major difficulty with the statistical evaluation of flares 
on different stars is that the sensitivity of the flare detection 
process depends on the measured (quiescent) count rate, which 
determines the signal-to-noise (S/N), and hence on the distance 
to the star The observational bias consists in the fact that for 
bright stars (Lqui large) the minimum luminosity Lp of a de- 
tectable flare is higher than for a faint star The result is, among 
others, that at first hand it can not be decided whether any ob- 
served correlation between Lqui and Lp is real or produced by 
this effect. In Fig. ^ we have plotted the flare luminosity Lp 
against the quiescent luminosity Lqui. 

The contribution of the observational bias to this correlation 
can be estimated as follows: For each quiescent count rate /qui 
we can determine the minimum strength Lp /Lqui needed for 
a flare to be detected, if we assume that a flare is found when- 
ever there is a rise in count rate of at least 3 a within one 400 s 
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is JD 24483 14. 125072 V47 1 Tau (200 1 07-0) 
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Time [h] 

= is JD 2448506.508770 LP357-4 (200556) 




r 



Time |h| 

-OisJD 2448498.598111 VBI41 (2007771 




Tims |hj 

= is JD 2448499.391016 VA275 (2007761 
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J___Ji 111, 1 



Timslhj 
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t - is JD 2448874.532536 VA677 (200945) 



^ IS JD 2448900.7881 13 VB16y(201097) 



-OisJD 2448875.989585 VA673 (201313) 
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Time |h| 

= IS JD 2448875.596344 VB14I(201368) 
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Fig. 4. PSPC lightcurves for flares on members of the Hyades cluster detected by the procedure described in Sect 4-. 1 . Identifi- 
cation of the X-ray source and 7?(95Ar observation request number (in brackets) are given for each source. Binsize is 400 s, 1 a 
uncertainties. The dashed Hne represents the quiescent count rate. Sohd Hnes are exponential fits to the data points belonging 
to the flare down to the quiescent emission. Background count rates are shown as upper limits when the background subtracted 
count rate is below zero. 
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Table 2. Properties of stars in Tau-Aur-Per which flared during at least one ROSATPSPC observation. Spectral types marked with 
an asterisk are determined from B — V or B — I given in the Open Cluster database using the conversion given by Schmidt-Kaler 
(1982). The meaning of the flags in column 'Multiplicity' is: B for binary, SBl for single-lined spectroscopic binary, and SB2 
for double-lined spectroscopic binary. For one flare on TTSs it is not clear whether it belongs to DD Tau, a cTTS, or CZTau, a 
wTTS, because the spatial resolution of the PSPC is too low to resolve these stars. DD Tau and CZTau both are binaries. Only 
flares on stars of spectral type G and later are analysed in the remainder of this paper. 



Designation 


Distance 


Sp.Type (Ref) 


vsini 


(Ref) 


Mult 


Binary sep 


TTS 






[pc] 






[km/s] 




["] 


Type 




T Tauri Stars 


LkHa270 


350.0 


K7M0 


(1) 










C 


tl) 


BPTau 


140.0 


K7 


(1) 


10.0 


(17) 






C 


(1) 


V410x-ray7 


140.0 


Ml 


(2) 










w 


12) 


HD283572 


140.0 


G51 


(1) 


75.6 


(18) 






w 


(1) 


DDTau/CZTau 


140.0 


Ml/Ml. 5 


(1) 






B/B 


0.57/0.33 (27) 


CAV 


(1) 


L1551-51 


140.0 


K7 


(1) 


27.0 


(1) 






W 


(1) 


RXJ0437. 5-1-1851 B 


140.0 


M0.5 


(3) 


10.5 


(3) 


B 


4.30 (28) 


W 


(3) 


LkCal9 


140.0 


KO 


(1) 


18.6 


(17) 






W 


(1) 


LHq92 


300.0 


KO 


(5) 










C 


tl) 


RXJ0422. 1-1- 1934 


140.0 


M4.5 


(4) 






B 


12.0 (28) 


C 


(4) 


TTau 


140.0 


KO 


(1) 


20.7 


(17) 


B 


0.71 (29) 


C 


(1) 


RXJ0255. 4-1-2005 


65.0 


K6 


(6) 


10.0 


(6) 






w 


(6) 


LkHa325 


350.0 


K7M0 


(1) 






B 


11.0 (30) 


w 


(1) 


Pleiads 


hiil384 


116.0 


A4V 


(7) 


215.0 


(19) 










hii2147 


116.0 


G9 


(8) 


6.9 


(20) 


SB2 


(31) 






hiillOO 


116.0 


K3 


(8) 


5.4 


(20) 


B 


0.78 (32) 






hii303 


116.0 


K1.4* 


(9) 


17.4 


(20) 


B 


1.81 (32) 






hii298 


116.0 


K1.4* 


(9) 


6.6 


(20) 


B 


5.69 (32) 






hcg307 


116.0 






13.0 


(21) 










hcgl44 


116.0 


















hii2244 


116.0 


K2.5 


(10) 


50.0 


(22) 










hcg422 


116.0 


















hiil516 


116.0 


K7* 


(11) 


105.0 


(20) 










hcgl43 


116.0 


M2* 


(11) 














sk702 


116.0 


















hiil74 


116.0 


K0.8* 


(9) 


28.0 


(23) 










hiil91 


116.0 


K7.7 


(12) 


9.1 


(20) 










hcg97 


116.0 


















hiil653 


116.0 


K6 


(8) 


21.0 


(8) 










hii345 


116.0 


G8 


(10) 


18.9 


(20) 










hcgl81 


116.0 


M1.5* 


(11) 














hii253 


116.0 


Gl 


(10) 


37.0 


(23) 










hii212 


116.0 


MO.O 


(12) 


10.0 


(8) 










Hyads 


V471Tau 


46.8 


K1.2* 


(13) 














VB50 


44.9 


GIV 


(14) 






SB? 


(33) 






VA677 


57.0 


K 


(15) 


24.0 


(24) 


SB2 


(24) 






LP357-4 


46.3 


M3 


(15) 














VA275 


46.3 


M2-3 


(15) 


10.0 


(25) 










VB141 


47.9 


A8* 


(16) 














VA334 


41.0 


MO 


(15) 


6.0 


(24) 










VB169 


46.5 


A6.6* 


(16) 














VA673 


46.3 


Ml 


(15) 






SB 


(34) 






VB190 


46.3 


K 


(15) 


8.5 


(24) 


SB 


(32) 






VB85 


41.2 


F5V 


(14) 


55.0 


(26) 




(35) 






VB45 


47.2 


Am 


(14) 


12.0 


(26) 


SBl 


(26) 







;VA-|vanAltena 196S|,LP-|Luyten etal. 198li hii -fiertzsprung 194^, hcg - ^aro et al. 198^, sk-S tauffer 
etal. 1991 

References: (1) - Herbig & Bell 1988, (2) - Strom & Strom 1994, (3) - Wichmann et al. (in preparation), (4) - Martin & Magazzii 1999, (5) - 
Herbig 1998, (6) - Hearty et al. 2000, (7) - Mendoza 1956, (8) - Stauffer & Hartmann 1987, (9) - Johnson & Mitchell 1958, (10) - Soderblom 
et al. 1993a, (11) - Stauffer 1982, Stauffer 1984, (12) - Prosser et al. 1991, (13) - priv. comm. between J. Stauffer and E. Weis according to 
the Open Cluster data base, (14) - Morgan & Hiltner 1965, (15) - Pesch 1968, (16) - Morel & Magnenat 1978, (17) - Hartmann & Stauffer 
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Table 3. Parameters derived from the lightcurves of flares on T Tauri Stars detected in the ROSAT PSPC observations from Table 1 (see Sect. 4.2 for an explanation). The 
last column gives the reference for flares that have been presented elsewhere in the hterature: [a] Strom & Strom (1994), [b] Preibisch et al. (1993), but newly reduced 
here. In the last two rows we give the mean (determined by taking account of lower/upper limits) and the median for each parameter. 



Desig. 


ROR 




-^max 




ec 


'^ris 


^qui 


-^'max 




Lp/Lqui 




Ef 


Notes 






[cps] 


[cps] 


[h] 


[h] 


[erg/s] 


[erg/s] 








[erg] 


[erg] 




LkHa270 


180185 


0.061 ± 0.008 


0.256 ± 0.008 


<8.85 




< 7.84 


8.13e+30 


> 3.41e+31 


> 


3.197 




2.59e+35 


5.23e+35: 




BPTau 


200001-0 


0.037 ± 0.004 


0.247 ± 0.034 


2.67 


+0.96 
-0.92 


> 0.22 


7.89e+29 


5.27e+30 




5.676 


± 1.110 


7.58e+33 


2.73e+34 




V410x-ray7 


200001-0 


0.050 ±0.011 


0.114 ±0.018 


0.80 


+7.22 
-0.53 


> 0.00 


1.07e+30 


> 2.43e+30 


> 


1.280 




3.07e+33 


> 2.52e+33: 




HD283572 


200001-1 


0.296 ± 0.069 


1.762 ±0.101 


1.78 


+0.18 
-0.19 


< 1.56 


6.31e+30 


3.76e+31 




4.953 


± 1.226 


4.04e+34 


1.27e+35 


[a] 


DDTau/CZTau 


200001-1 


0.004 ± 0.005 


0.113 ±0.013 
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Table 4. Parameters derived from the lightcurves of flares detected on members of the Pleiades cluster in the ROSAT VSPC observations from Table 1. The last column 
gives the reference for flares that have been presented elsewhere in the literature: [a] Gagne et al. (1995), but newly reduced here. In the last two rows we give the mean 
(determined by taking account of lower/upper limits) and the median for each parameter. 
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Table 5. Parameters derived from the lightcurves of flares detected on members of the Hyades cluster in the ROS AT PSPC observations from Table 1 . No analysis of flares 
on Hyads are reported in the literature. In the last two rows we give the mean (determined by taking account of lower/upper limits) and the median for each parameter. 
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Fig. 5. Correlation between quiescent and flare luminosity. 
TTSs, Pleiads, and Hyads are represented by different plotting 
symbols: TTS - circles, Pleiads - triangles, Hyads - crosses. 
The solid curve is the minimum flare strength needed for de- 
tection of a flare with given Lqui if the star is at a distance of 
140 pc. The data points below that curve represent no contra- 
diction to the calculated threshold because the detection thresh- 
old increases with stellar distance, i.e. apparent brightness of 
the stars, and therefore Pleiads and Hyads may show flares with 
smaller ip for given Lq^i- 



time bin. (In our actual flare search we were even more con- 
servative; see Sect. ^.) Hypothetical events of that kind obey 
a detection threshold curve for Lp/iqui as shown in Fig. ^. 
As mentioned above, the minimum flare luminosity needed for 
detection of a flare becomes larger with increasing quiescent 
brightness. In contrast, the required luminosity ratio, i.e. the 
relative strength of the events, decreases when /qui increases. 
Note also, that the curve in Fig. ^ is distance independent. But 
the relation between Lqui and the corresponding minimum Lp 
of a detectable flare differs for stars at different distances. In 
Fig. H we have overplotted the theoretical threshold for detec- 
tion of a flare on a star at 140 pc distance. Note, that the slope 
of the data in Fig. H is somewhat steeper than the increase of 
the threshold imposed by the S/N. This seems to indicate an in- 
trinsic correlation between quiescent and flare luminosity. We 
have subtracted the theoretical threshold value for Lp from the 
observed flare luminosity for each of the stars from Fig. |[ Cor- 
relation tests for the difference between threshold and observed 
value for Lp, (Lp.theo — -Z^F.obs), with Lqui show that the cor- 
relation is of low significance, q;=0.05. The data points below 
the theoretical curve are all Pleiads or Hyads. They do not con- 
tradict the threshold curve, since Pleiades and Hyades stars are 
closer than 140 pc and therefore have a lower flare detection 
threshold. 




-3 -2 -1 

Ig (Quiescent Count Rate I^^j) [cps] 

Fig. 6. Relative flare strength, Lp/Lqui, as a function of the 
quiescent count rate /qui in double logarithmic scale. Assum- 
ing that the flare is characterized by a rise in count rate of 3 cr 
above the mean during a single 400 s time bin, the solid curve 
gives the minimum values of Lp / Lqui for which a flare will be 
detected on stars with given /qui. The observed values consis- 
tently lie above this line. The meaning of the plotting symbols 
is the same as in Fig. |[ For clarity we have omitted error bars. 



6. Statistical comparison of the flaring stars 

We present now a statistical analysis of the X-ray flares from 
Tables ^ ^ and |[ A detailed discussion of the (quiescent) X- 
ray properties of all detected and undetected stars is postponed 
to a later paper (Stelzer et al., in preparation). 

In this section, different flare parameters will be checked 
for dependence on age, circumstellar environment, and rotation 
rate (Sect. x2,^J, ^.) to see whether any of these properties 
has an effect on the characteristic luminosity and time scales 
of coronal activity. For the statistical comparison of the flar- 
ing stars the ASURV package version 1.2 (Feigelson & Nelson 
1985) was used. 

First we compare the flaring populations of TTSs, Pleiads, 
and Hyads concerning their effective temperatures. We have 
converted spectral types to effective temperatu res using the 
conversion given in Kenyon & Hartmann (1995) for PMS stars 
earlier than MO, and t^uhman (1999) for PMS M-type stars in- 
termediate between dwarfs and gian ts. For Pleiades and Hya des 
stars we have used the conversion of Schmidt-Kaler (1982) . We 



have applied two-sample tests to each pair of Tcff-distributions 
to reveal possible differences between flaring stars of the three 
groups. Henceforth, we denote the probability that the distri- 
butions are similar by a. In all but one of the comparisons 
we found a > 0.2, and therefore no significant differences in 
Tcff . The exception is the logrank test between TTSs and Hyads 
where a — 0.03. 

Most flares occurred on G, K and M stars. However, some 
events were observed on A and F stars. Stars of intermediate 
spectral type, lacking both a convection-driven dynamo and 
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Table 6. Results of two-sample tests between each pair of stel- 
lar samples (TTSs, Pleiads, and Hyads) with respect to the ef- 
fective temperature Teg of the flaring star. Only flares on G, 
K, or M stars have been admitted. The analysis was performed 
with the ASURV package. Next to the mean and median of 
Toff we give the probability that the null hypothesis of two dis- 
tributions being the same is true derived from Gehan's general- 
ized Wilcoxon test and the logrank test. The values suggest that 
there is no significant difference between the spectral types of 
the flaring TTSs, Pleiads, and Hyads and justify to combine G, 
K, and M stars for the comparison of flares from these different 
stellar groups. 

Ig Toff- 



Sample 
Name 


Sample 
Size 


MEAN 


MED. 


Prob. 
GW (HV) 


Prob. 

logrank 


TTSs 
Pleiads 


15 
14 


3.63 ± 0.08 
3.66 ± 0.07 


3.600 
3.675 


0.395 


0.525 


TTSs 
Hyads 


15 
11 


(see above) 
3.63 ± 0.07 


3.595 


0.897 


0.972 


Pleiads 
Hyads 


14 
11 


(see above) 
(see above) 




0.212 


0.401 



a strong stellar wind, seem to have no efficient mechanism 
to generate X-ray flares. Therefore, it is often assumed that 
X-ray emission apparently seen on A or B stars, can be at- 
tributed to an (unknown) late-type companion (see e.g. Stauf- 
fer et al. 1994, [Gagnc & Caillault 1994| , [Panzera et al. 1999| ). 
The same arguments can be applied to explain X-ray flares on 
these stars. In any case, the emission mechanism of early-type 
stars is different from that of late-type stars. From the sharp 



onset of rotation-activity relations in dwarf stars Walter (1983) 
has argued that the onset of solar-like dynamo activity occurs 
abruptly at about spectral type F5. To ensure that no stars with 
X-ray generation mechanisms other than stellar dynamos are 
included, we have excluded the stars of spectral type F and ear- 
lier from the statistical analysis presented in this paper, i.e. we 
have restricted the flare sample to events on G, K, and M stars. 
This limitation provides samples which have similar Tcff distri- 
butions, i.e. a > 0.2 also for the two-sample test between TTSs 
and Hyads (see Table ^ for the detailed results). This justifies 
to combine all flaring late-type stars for the statistical analysis. 
In the following the stellar sample is restricted to G, K, and 
M stars. If not explicitly mentioned the two flares on known 
white-dwarf systems (on V471 Tau and VA673) are excluded 
from the sample, since the white dwarf could be responsible for 
the X-ray event instead of its late-type companion. 

6.1. Flare frequency of MS stars and spectral type 

It is interesting to ask whether the depth of the convection 
zone has any influence on the occurrence of surface flares. 



Since the relative size of the convection zone increases for later 
spectral types, the distribution of flares onto stars of different 
spectral types may help to solve this question. What we really 
want to check is whether the flare frequency depends on stel- 
lar mass, which corresponds to spectral type on the MS. Be- 
cause PMS stars still evolve through the Hertzsprung-Russell 
diagram (HRD), i.e. change their spectral type, we exclude 
the TTSs from this part of the analysis. Flaring Pleiades and 
Hyades stars are combined to increase the sample size. 

We have studied the spectral type distribution of flares by 
comparing the number of flares on stars of a certain spectral 
type to the total number of detected stars of that spectral type. 
The detection sensitivity for flares of a given strength L-p / Lqui 
is different for each star because it depends on the level of qui- 
escent emission /qui (see Sect. ||). /qui, depends on the spec- 
tral type of the star For this reason, a simple comparison be- 
tween numbers of flares and numbers of detected stars of each 
spectral type would be misleading. The observational bias can, 
however, be eliminated if the numbers (of flares and detections) 
are evaluated above a certain threshold Lp/iqui- We compare 
the number of flares with measured luminosity ratio above a 



critical value (Lp/L, 



qui^cnt 



to the number of detected stars 



for which /qui exceeds the minimum value needed for detec- 
tion of a flare of that critical strength. We have compiled these 
numbers for a reasonable range of values Lp/iqui, and show 
the result in Fig. 0. Plotted are the number of flares exceeding 
Lf I -^qui divided by the number of detected stars that are bright 
enough for detection of flares with that value of Lp/Lqui. G 
stars clearly show the smallest rate of events throughout all of 
the observed range of flare strengths. 



6.2. Age of flaring stars {Luminosity functions) 

To study how the flare activity of young late-type stars evolves 
with stellar age we have computed luminosity distribution 
functions (LDF) and performed two-sample tests for three sub- 
samples of stars: TTSs, Pleiads, and Hyads. 

Maximum likelihood distributions for TTSs, Pleiads, and 
Hyads are presented in Fig. ^ for both flare luminosity Lp and 
mean luminosity during the quiescent part of flare observations 
Lqui. Note, that Fig. ^ (b) contains no upper limits because 
only stars which have shown a flare are included, and iqui dur- 
ing flare observations can be extracted from the lightcurves. 
The flare luminosity in Fig. ||(a) includes upper limits. Since 
ip Lmax - Lqui, Upper Umits for L,„ax (see Table |, |, 
and U) translate to upper limits for Lp. LDFs for all non-flaring 
stars (detections and non-detections) will be shown elsewhere 
(Stelzer et al., in preparation). 

Two-sample tests were applied to each pair of LDFs to 
search for differences. The results are given in Table | (for 
Lp and Lqui). The null hypothesis of two samples being the 
same is rejected for all pairs of flare luminosity distributions 
at significance levels a < 0.05. The quiescent luminosity of 
flaring TTSs is different from both the quiescent luminosity of 
flaring Pleiads and flaring Hyads. Usually the quiescent lumi- 
nosity functions of Pleiades and Hyades stars are also found to 
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Fig. 7. Fraction of flares among detected stars above a given 
threshold Lp/Lq^i for different spectral types. Shown are G 
stars (dashed line), K stars (solid line) and M stars (dotted 
line). Pleiades and Hyades stars have been combined. TTSs are 
not considered since they are on the PMS where the relation 
between the stellar interior structure and the X-ray emission 
might be different. Throughout the examined range of Lp / iqui 
the number of flares per detected stars is lowest for spectral 
type G. 



Table 7. Results of two-sample tests for the flare luminosity Lp 
and quiescent luminosity of flare observations Lqui- The num- 
ber of upper limits to the luminosity are listed in brackets in 
column 'Sample Size'. See Table ^ for the meaning of the re- 
maining columns. All distributions of Lp are significantly dif- 
ferent from each other. Flaring Pleiads and Hyads have similar 
Lqui distributions. 

Sample Names Sample sizes Prob. GW Prob logranlc 



Flare luminosity Lp 



TTSs - Pleiads 
TTSs - Hyads 
Pleiads - Hyads 



15(6) 
15(6) 
14(5) 



14(5) 
9(3) 
9(3) 



0.021 
0.001 
0.007 



0.017 
0.000 
0.010 



Quiescent luminosity of flaring stars Lqu 



TTSs - Pleiads 
TTSs - Hyads 
Pleiads - Hyads 



15(0) 
15(0) 
14(0) 



14(0) 
9(0) 
9(0) 



0.025 
0.007 
0.615 



0.010 
0.001 
0.243 



be distinct (see e.g. Zaillault 1996). However, we find no dif- 



ference (a > 0.61) between the quiescence luminosities of the 
flaring stars in these two clusters. Using the relation between 
/qui and the threshold for Lp/Lq^i (see Fig. ^ we have de- 
termined that more than 90% of the detected Hyades stars are 
bright enough for detection of a flare whose strength Lp / Lq^i 
is equal to the mean observed for flares on late-type Hyads, i.e. 
Lp/Lqui = 4.527. The fact that mostly X-ray bright Hyades 
stars display flaring activity is therefore not a selection effect. 
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□ Hyades 
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29 30 31 
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Fig. 8. Luminosity functions for the flaring TTSs, Pleiads, and 
Hyads during the X-ray flare (a) and during quiescence (b). 
Stars of spectral type earlier than G or unknown spectral type, 
and known white-dwarf systems are excluded. Note that the 
quiescent distribution of the Hyades stars is consistent with that 
of the Pleiads, due to the selection of flaring stars only (see 
discussion in the text). 



Instead, flaring Hyads indeed are overluminous compared to 
the non-flaring Hyades stars detected by the ROSAT PSPC. 

The mean luminosities of TTSs, Pleiads, and Hyads and 
their standard deviations derived with inclusion of upper limits 
are given in Table 

6.3. Flaring cTTSs and wTTSs 

So far we have not distinguished between cTTSs and wTTSs, 
because it is a matter of debate whether all cTTSs are younger 
than wTTSs. However, they are clearly distinguished by their 
circumstellar environment. The disks of cTTSs may influence 
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Table 8. Mean and median flare and quiescent luminosity of 
flaring TTSs, Pleiads, and Hyads. The sample of cTTSs was to 
small to compute the median of Lp. 



Sample Name 


Sample Size 


Mean 


Median 


Flare Luminosity Ig Lp 


TTSs 


15 


(6) 


31.05 ±0.19 


30.65 


Pleiads 


14 


(5) 


30.51 ±0.12 


30.30 


Hyads 


9 


(3) 


30.06 ±0.11 


29.99 


cTTSs 


6 


(4) 


31.44 ±0.32 




wTTSs 


8 


(1) 


30.81 ±0.22 


30.58 


Quiescent Luminosity of Flaring Stars Ig 


^qui 


TTSs 


15 


(0) 


29.98 ±0.17 


29.92 


Pleiads 


14 


(0) 


29.52 ±0.11 


29.61 


Hyads 


9 


(0) 


29.50 ±0.07 


29.47 


cTTSs 


6 


(0) 


29.94 ±0.19 


29.82 


wTTSs 


8 


(0) 


30.22 ±0.19 


30.04 



Table 9. Results of two-sample tests for differences between 
cTTSs and wTTSs with respect to the flare parameters Lqui, 
Lp, Tdoc, and Lp/iqui- The flare on the unresolved stars 
DDTau/CZTau is not considered here, since it is not clear 
whether the event should be attributed to the wTTS (CZ Tau) 
or to the cTTS (DD Tau). The meaning of columns 2-5 is the 
same as m Table 



Parameter 


Sample Size 


Prob. GW 


Prob logrank 




cTTSs 


wTTSs 






Ig Lp 


6(4) 


8(1) 


0.053 


0.046 


Ig iqui 


6(0) 


8(0) 


0.206 


0.264 


"^dcc 


6(2) 


8(1) 


0.714 


0.752 


Lf / Lqui 


6(4) 


8(1) 


0.007 


0.002 



flare activity. We have, therefore, compared cTTSs and wTTSs 
with respect to several flare parameters (see Table || for the 
results). Significant differences are found in the flare luminos- 
ity Lp and relative strength of the flare Lp/^qui- The decay 
timescale does not depend on the type of TTS. The values given 
in Table ^ have been derived from all flares on TTSs except the 
one on DD Tau / CZ Tau. DD Tau is a cTTS binary and CZ Tau 
a wTTS binary. The four stars are not resolved in the PSPC 
image. It is therefore impossible to classify this flare concern- 
ing the type of TTS. We have performed two further series of 
two-sample tests in which the event is included. In one of these 
series of tests the flare is attributed to DD Tau, and the other 
time to CZ Tau. The significance of the results did not change. 

The mean flare and quiescent luminosities of cTTSs and 
wTTSs are given in Table |^. cTTSs, although characterized by 
lower quiescent emission, show stronger flares than wTTSs. 



6.4. V sin i of the flaring stars 

Stellar rotation is one of the necessary conditions for magnetic 
activity. We have, therefore, examined the influence of the stel- 
lar rotation rate on the characteristics of X-ray flares. The rela- 
tion between flare parameters and projected rotational velocity, 
V sini, is shown in Fig. 0. The statistical significance for cor- 
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Fig. 9. X-ray luminosities versus projected rotational velocity 
V smi in double logarithmic scale: top - flare luminosity Lp, 
bottom - luminosity ratio Lp/Lqui- While Lp and v s'mi may 
show a weak positive correlation, Lp / Lqui does not depend on 
the rotation. 



relations between some flare parameters and v s'mi is given 
in Table |o| (columns 2 and 3). The weak correlation between 
luminosity, both Lp and iqui, and v s'mi is significant. The 
decay time Tdoc and the relative flare strength Lp / iqui, on the 
other hand, are not related to v sin i. 

We have studied the flaring population in terms of differ- 
ences in flare characteristics between slow and fast rotators. 
The boundary was set to 20 km/s because this choice gives 
two samples of about equal size: 16 slow and 12 fast rota- 
tors showed an X-ray flare. Stars from Tables ^ ^ and ^ for 
which no measurement of v sin i is available are ignored. The 
result of two sample tests for the parameters Lp, Lqui, T'doc, and 
Lp /Lqui are presented in the remaining columns of Table |l^ 
In no case the null hypothesis that slow and fast rotators are 
drawn from the same distribution was rejected at significance 
level a < 0.05. 
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Table 10. Relation between flare parameters Lqui, Lp, Tdcc, 
and Lp/Lqui and the projected rotational velocity v s'mi: 
Columns 2 and 3 are results of correlation tests between each of 
the listed parameters and v sini. The remaining columns give 
results of two-sample tests for differences between slow and 
fast rotators with boundary at 20 km/s. The samples consist of 
16 slow versus 12 fast rotators. 



Parameter 


Pearson r 


Spearman r 


Prob. 


Prob 




(Signif.) 


(Signif.) 


GW 


logrank 


Ig Lp 


0.479 (0.005) 


0.372 (0.026) 


0.075 


0.091 


Ig iqui 


0.468 (0.006) 


0.496 (0.004) 


0.269 


0.528 


^dcc 


0.247 (0.103) 


0.221 (0.129) 


0.487 


0.584 


Lf / Lqui 


-0.017 (0.466) 


-0.287 (0.069) 


0.399 


0.547 



7. Flare rates 

In this section, we will derive flare rates as a means to deter- 
mine the activity level for a stellar sample with distinct proper- 
ties. The characteristic properties which will be examined are 
(a) stellar age (comparing TTSs, Pleiads, and Hyads), (b) stellar 
rotation (comparing slow and fast rotators), and (c) stellar mul- 
tiplicity (comparing close binaries to other stars). Flare rates 
will be computed separately for each group of stars. 

We assume that the duration of the active state is repre- 
sented by the decay timescale Tdcc, i e- the generally poorly 
restricted rise times Tiis are neglected. This is certainly wrong 
for the flare on hcg 144 which seems to have a somehow re- 
versed character (slow rise and rapid decay). However, hcg 144 
is a star of unknown spectral type and therefore not part of 
the group to be examined here. To compile the flare rates, 
F = {Ti)/Tc,hs, we have added up the decay timescales Tdoc 
of the flares and divided this sum by the total observing time, 
Tohs, of all detections (flaring and non-flaring stars). Only the 
nearest identification of each X-ray source has been considered 
in the compilation of Tobs- But for multiple systems we have 
multiplied the observing time by the number of components. 
For the compilation of Tobs we have eliminated data gaps larger 
than 1 h, the typical flare duration. This provides us the fraction 
of the total observing time during which the stars are observed 
in the active state. 

In practice ^ (r^) is computed from the sample mean f 
returned by ASURV's Kaplan-Meier estimator. This way we 
ensure that upper limits to Tdcc are taken into account. The 
Kaplan-Meier estimator returns also the uncertainty of f. To 
include the spread of the data in the estimation of F we have 
converted this uncertainty of the mean to the sample variance 
ar- Consequently: 

^^f_N^._^ (1) 

Tobs Tobs 

7.7. Flare rate and stellar age 

The evolution of flare rates with stellar age is examined by 
comparing the flare frequency of TTSs to that of the Pleiades 
and the Hyades. 15 flares have occurred on TTSs, 14 on late- 



type Pleiads, and 1 1 on late-type Hyads. We have derived the 
foflowing values for the flare rate F: 0.86 ± 0.16% (TTSs), 
0.67 ± 0.13% (Pleiads), and 0.86 ± 0.32% (Hyads). When the 
white-dwarf binaries are excluded the flare rate for Hyades de- 
clines to Tn^noWD = 0.71 ± 0.30%. 

The flare rates are biased for several reasons which will 
be explained next. First, the flare detection limit is determined 
by the S/N, which in turn depends on the distance to the star. 
Therefore, flare rates of TTSs, Pleiads, and Hyads are only 
comparable above a limiting minimum strength of the flare, 
expressed by a threshold Lp/Lqui- And, secondly, incomplete 
data sampling might lead to wrong conclusions about the decay 
timescale of individual events and thus contaminate the result- 
ing F. 

To solve the first problem we have scaled the quiescent 
count rate of all flaring stars to a distance of 140pc, the dis- 
tance of most of the TTSs. I.e. we have multiplied all quiescent 
count rates with a factor (r-r^ — These theoretical values of 

V 14u pc ' 

/qui correspond to higher flare detection thresholds Tp / Lqui 
for all stars except the ones in Perseus. All flares from Perseus 
stars would be detected at 140pc, since they are further away 
than this distance. The observed luminosity ratios of all flar- 
ing stars have then been compared to the theoretical thresh- 
old needed if the star were at 140 pc. All flares for which the 
observed value is below this requirement should be neglected 
when the flare rates are computed. It turns out that all flares on 
Pleiads remain above the 140 pc threshold. But only 7 out of 
1 1 flares on Hyads (one of the 7 is a white-dwarf system) have 
Tp/Tqui high enough to be detected at a distance of 140 pc. 
Now the comparison of our different samples is free from the 
sensitivity bias. And we derive flare rates F of 0.86 ± 0.16% 
(TTSs), 0.67 ± 0.13% (Pleiads), and 0.46 ± 0.19% (Hyads). 
Without the white-dwarf binary F decreases for the Hyades to 
0.32 ±0.17%. 

The uncertainties in the measurement of the flare duration 
are less easy to overcome. The large flare rate of TTSs is par- 
tially due to two extraordinary long events of duration > 8 h 
(see Table The decay times of both of these flares are con- 
sidered to be an upper limit. If these two events are discarded 
from the sample of flares, Ttts = 0.74 ± 0.14%. 

We have also compiled F for cTTSs and wTTSs separately 
to see whether the circumstellar environment has any influence 
on the frequency of the flare activity. Among the events on 
TTSs, 6 are observed on cTTSs and 8 on wTTSs. An additional 
flare was seen from the unresolved stars DDTau/CZTau. 
The classification of this event within the subgroups of TTSs 
remains therefore unclear, and complicates the comparison 
of F for the two classes of TTSs. At first, the event on 
DDTau/CZTau has been eliminated from the sample, thus 
that 6 flares on cTTSs are opposed to 8 flares on wTTSs. 
The respective flare rates are F^ = 1.09 ± 0.39% and F^ = 
0.65±0.16%. When the ambiguous event is counted on the side 
of the cTTSs F^ rises to 1.28 ± 0.37%. When it is attributed to 
the wTTS CZTau instead, F^ becomes 0.76 ± 0.16%. Note, 
that even though the number of flares on wTTSs is higher than 
the number of flares on cTTSs, the flare rate for wTTSs is lower 



B. Stelzer, R. Neuhauser & V. Hambaryan: X-ray flares in Taurus-Auriga-Perseus 



19 



<> 

-<>- 



1.4 
1.2 
1.0 

^ 0.6 

0.4 

0.2 
0.0 

5 6 7 8 9 

log Age [yrs] 

Fig. 10. Flare rate F = {Ti)/Tohs as a function of stellar 
age for TTSs, Pleiads, and Hyads. To eliminate biases due to 
the S/N-dependence of the detection sensitivity for flares, only 
events which are bright enough for detection at a uniform dis- 
tance of 140 pc are considered. Furthermore, flares which can 
not definitely be assigned to a late-type star are excluded. The 
flare rates for cTTSs (diamonds) and wTTSs (triangles) alone 
are overplotted. The upper of the two symbols applies if the 
ambiguous event on DDTau/CZTau is included in the sam- 
ple (see text). The horizontal bars represent the age spread of 
TTSs. 



than the flare rate for cTTSs. This is possible because of differ- 
ences in the total observing time. 

F as a function of stellar age is displayed in Fig. [l^. The 
decline of the flare rate with stellar age is obvious. Rates for 
cTTSs and wTTSs are symbolized by diamonds and triangles, 
respectively. The location of the lower diamond and triangle 
describes the flare rates without the event on DDTau/CZTau. 
The upper diamond and triangle are values for F if this flare is 
included in the respective group of TTSs. 

7.2. Flare rate and rotational velocity 

Here we examine whether the flare frequency depends on ro- 
tation. This is done by computing F (defined as before) for 
fast rotators on the one hand {v sini > 20km/s) and slow 
rotators on the other hand {v sini < 20km/s). Again, only 
late-type stars are considered. The resulting rates are Fgiow = 
0.55 ± 0.10% and Ffast = 1-55 ± 0.38%. Thus there is a clear 
trend towards an increase of flare activity with increasing rota- 
tional velocity. 

7.3. Flare rate and multiplicity 

Another interesting question is whether the circumstellar sur- 
roundings have any influence on the flare frequency. The coro- 
nal activity may e.g. change if there are interactions between 



the magnetic fields of binaries. Such interactions are expected 
to take place only in close binaries. To search for such a con- 
nection we, therefore, discriminate between spectroscopic bi- 
naries on the one hand and all others, i.e. singles or visual mul- 
tiples. The flare rate F is computed in the same way as be- 
fore. Since the observation time of each stellar system has been 
multiplied by the number of components, the flare rates should 
be about equal for both samples if the underlying physics 
are the same. However, we find that the flare rate of spec- 
troscopic binaries is enhanced by more than a factor of two: 
Fnon-SB = 0.64 ± 0. 12% and FgB = 1.43 ±0.25%. Note, fliat 
the study of individual flare parameters (similar to the analysis 
of Sect. g|) has shown no difference for these two samples. 

8. Hardness Ratios 

For most of the flaring sources not enough counts are collected 
by the PSPC to compare the different levels of X-ray emission 
in a detailed spectral analysis. Therefore, we use hardness ra- 
tios to mark spectral changes. ROSAT PSY'C hardness ratios are 
defined by: 



HRl = 



H_S 

H + S 



HR2 = 



H-2-H, 



(2) 



where S, H, Hi, and H2 denote the count rates in the 
ROSAT PSPC soft (0.1-0.4keV), hard (0.5-2.0keV), hardl 
(0.5-0.9 keV) and hard2 (0.9-2.0keV) band respectively. For 
each flare observation HR 1 and HR 2 are computed for three 
activity stages representing the quiescent state (pre- and post- 
flare), the rise and the decay, respectively. Sometimes no counts 
are measured in one or more of the energy bands. Whenever 
this is the case we have derived upper limits for the hardness 
ratio making use of the background counts in that energy band 
at the source location. 

The observed hardness ratios, HR 1 and HR 2, are plotted 
in Fig. [n]. The plots comparing quiescent and flare state show 
marginal evidence that most of the stars lie below the diagonal 
in the hardness plot (see lower left panel of Fig. |ll]) and thus 
are harder during the flare intervals as compared to their qui- 
escence. No significant difference in hardness is observed be- 
tween flare rise and flare decay. When impulsive heating takes 
place before the outburst the plasma cools quickly by radia- 
tion and conduction to the chromosphere. Therefore, the simi- 
lar hardness observed during rise and decay phase suggests that 
heating takes place throughout the decay. 

To quantify the differences in hardness between different 
flare stages we have computed mean hardness ratios for each of 
the stellar groups. In Table [u] we show the mean hardness for 
each activity stage (quiescence, rise, and decay) and each sam- 
ple of stars (TTSs, Pleiads, and Hyads). The hardness changes 
systematically when the three groups are compared to each 
other: TTSs display the hardest spectra, followed by Pleiads, 
which in turn are characterized by higher hardness ratios than 
the Hyades stars. This is also manifest in the hardness plots of 
Fig. Owhere the three samples occupy different regions. In 
Sect. |6.2| it was shown that the flare luminosity declines with 
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Fig. 11. Comparison of ROSAT PSPC hardness ratios for quiescent state, flare rise and flare decay of the late-type stars from 
Tables |, |, and |} top - HR 1, bottom - HR 2. The meaning of the plotting symbols is the same as in Fig. ^. Stars located on the 
dashed line show no change in hardness between quiescent and flare state. 



stellar age. As a consequence, the spectral hardness and the 
flare luminosity are correlated. The relation between hardness 
ratios and Lp is displayed in Fig. |l2|and suggests that the more 
luminous flares are associated with hotter plasma. 

9. Discussion 

9.1. Methods for flare detection 

Using binned data to detect flares introduces observational re- 
strictions. The sensitivity for detection of small flares is lower 
and very short flares remain unobserved due to the time bin- 
ning. Apart from these limitations, our flare detection produces 
reliable results as verified by comparison to both an alternative 
approach using Bayesian statistics and, where possible, previ- 
ous detections of flares by visual inspection stated in the litera- 
ture. 

The importance of Bayesian statistics to astronomical time 



series analysis has been described by Scargle (1998) and first 
appUed to ROSAT observations of flare stars by Hambaryan 
et al. (1999). This approach, unlike the 'classical' method used 
here, works on the raw, unbinned data and therefore has a time 
resolution which is only limited by the instrument clock. 

We have performed a detailed comparison of the events rec- 
ognized by the two methods. For the flare detection with the 
Bayesian algorithm the prior odds ratio, Opii, was set to 1 . This 







HRl 






Quiescence 


Rise 


Decay 


ITS 


0.75 ±0.38 


0.54 ±0.47 


0.70 ± 0.44 


Pleiads 


0.54 ±0.26 


0.50 ±0.29 


0.58 ±0.30 


Hyads 


-0.09 ±0.24 


-0.16 ±0.15 


-0.07 ±0.16 






HR2 






Quiescence 


Rise 


Decay 


ITS 


0.36 ±0.24 


0.38 ±0.31 


0.39 ±0.21 


Pleiads 


0.07 ±0.20 


0.08 ±0.22 


0.21 ±0.14 


Hyads 


-0.01 ±0.12 


-0.04 ±0.27 


0.14 ±0.27 



Table 11. Evolution of the ROSATPSPC hai'dness ratios HR 1 
and HR 2 during stellar flares. Given are the weighted means 
of flares on TTSs, Pleiads, and Hyads measured during three 
different phases selected from the X-ray lightcurve: Quies- 
cence (pre- and post-flare). Rise, and Decay. A slight increase 
in hardness is observed between quiescence and rise phase, but 
is not significant due to the large spread of the data within one 
sample. During the decay the spectrum retains its hardness. 
This indicates additional heating taking place during the de- 
cay phase. For each of the flare stages the spectral hardness 
decreases from TTSs over Pleiads to Hyads. 



means that at the beginning one-rate and two-rate Poisson pro- 
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Fig. 12. Correlation between ROSATVSiPC hardness ratios and 
flare luminosity L-p. The meaning of the plotting symbols is 
the same as in Fig. 1 1 Lower limits for Lp are shown when 
a data gap proceeds the observed maximum of the lightcurve 
(see Tables HI and I) 



cesses are assumed to have the same probability for being the 
correct description of the data set. The significance of any de- 
tection of variability is then given by the value of the posterior 
odds ratio, O21 ■ Applied to our data, 62 events are found at the 
5 fj level, and 95 events have O21 corresponding to at least 3 a. 
All but 5 of the flares discussed in this paper were among the 
5 fj detections. The remaining ones are detected at > 3 cr. But 
note that with the Bayesian method we find variabiUty in 182 
lightcurves (in contrast to our 52 flares). 

Although the Bayesian approach is sensitive to short events. 



we have persisted on the criteria explained in Sect. 4. 1 for two 



reasons: (i) While Bayesian statistics are sensitive to all kinds 
of temporal variability, we are here interested in large flare 
events only. This makes an additional selection process nec- 
essary, (ii) Comparison with the classical flare search used in 
this paper has shown that the Bayesian method needs further 
refinement. E.g. the outcome of the present algorithm used to 
search for flares depends sensitively on the value of the prior 
odds ratio. 



9.2. Interpretation of the results 

Before flares on different stellar groups are compared, it must 
be checked whether the composition of these samples is simi- 
lar The X-ray luminosity of MS stars depends on their spectral 
type. Therefore it would be desirable to separately investigate 
the flare activity from stars of different spectral types. However 
this is hindered by the low flare statistics. We have performed 
statistical tests where the flaring TTSs, Pleiads, and Hyads have 
been compared regarding to their Toff and hence spectral types. 
These tests have shown that it is justified to jointly analyse 
flares on all late-type stars, i.e. stars with spectral type G, K, 
andM. 

We have shown that the relative number of flares increases 
when going from spectral types G to K (see Fig. |^). Hereby, we 
have taken into account that the detection sensitivity for flares 
depends on the level of measured quiescent emission and hence 
on the spectral type. An interpretation is that deeper convection 
zones are favorable to the occurrence of surface flares. 



9.2.1. Age 

We found that in terms of absolute flare luminosity and en- 
ergy output TTSs surpass both Pleiads and Hyads. The mean 
flare luminosity of TTSs (Lf.tts = 1-13 10'^^ erg/s) is almost 
an order of magnitude higher than that for Hyads (Lp^Hya ~ 
1.15 lO'^" erg/s). The mean Pleiades flare luminosity is inter- 
mediate between that for TTSs and Hyades stars with ip^pio = 
3.26 10"^" erg/s (see also Fig. ^ and Table|^). This is partly due 
to the different distances of our stellar samples which result in 
different detection sensitivities for flares. Note, however, that 
this effect can explain only why no events with small L-p are 
observed on TTSs. But the lack of large events on Hyades stars 
is real. In Sect. ^ flare rates for TTSs, Pleiads, and Hyads have 
been established from an evaluation of the observed flare du- 
rations and the total observing time. Both, flare rate and mean 
flare luminosity decline with increasing stellar age. 

The quiescent luminosity of Hyades stars which showed a 
flare is larger than the average Lqui of Hyads (see Sect. 5.2). 
More than 90% of the detected Hyades stars are bright enough 
for detection of an average Hyads flare. Therefore, this result is 
not a selection effect, and we can conclude that only the most 
X-ray luminous Hyades stars exhibit X-ray flares. The interest- 
ing question whether the enhanced X-ray luminosity of flaring 
Hyades stars can be explained by their rotation rate can not be 
pursued with this set of data, because only for half of the flaring 
Hyades stars measurements of v sin i are available. 



9.2.2. Circumstellar Envelope 

If magnetic interactions between star and disk take place, the 
field lines will constantly b ecome twisted by differential rota- 
tion ( Montmerle et al. 2000 ). This may provide an environment 
favorable for magnetic reconnection and related flare activity. 

Six of the observed flare events can be attributed to cTTSs 
and 8 events to wTTSs. One of the flares on TTSs occurred ei- 
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ther on DD Tau, a cTTS, or on CZ Tau, a wTTS, both of which 
are not resolved in the ROSATPSPC observations. Two-sample 
tests show clear indications that flares on cTTSs are more X- 
ray luminous than those on wTTSs (see Table This holds no 
matter on which side the ambiguous event is counted. The flare 
rate is also slightly higher for cTTSs than for wTTSs, how- 
ever with low significance. Given the fact that quiescent X-ray 
emission of wTTSs is stronger than in cTTSs, this observation 
is surprising. A possible interpretation is that the stronger flare 
events on cTTSs may be due to violent interaction with their 
disks. 

9.2.3. Multiple Flares 

During four observations a second flare followed the first one 
(see lightcurves of VA334, VB 141, RXJ 0437.5H-1851B, and 
TTau in Figs. ^ and 0). From the number of observed flares 
and the total observing time the average duration between two 
flare events is estimated to be > 100 h. Therefore, from a sta- 
tistical point of view it is very unlikely to observe so many un- 
related 'double events'. We note, that double flares have been 



reported in the optical. And puenther & Ball (1999)| have pre- 
sented two flares that occurred within a few hours from the 
wTTS V819Tau. 

A possible interpretation of multiple flares is the star-disk 



scenario proposed by Montmerle et al. (2000) and mentioned 



in the previous subsection. However, this model does not seem 
to be accurate for our objects, which are more evolved and in 
part are known not to possess disks. 

9.2.4. Projected Rotational Velocity 

The statistical tests we have performed to discriminate between 
slow and fast rotators (with boundary drawn at 20km/s) reveal 
no dependence of individual flare parameters Lp, Lqui, Tdoc, 
and Lp / Lqui on the rotation rate. However, the flare frequency 
is about three times higher for fast rotators as compared to slow 
rotators: F^iow = 0.55 ± 0.10% and Ffast = 1-55 ± 0.38%. 

9.2.5. Binary Interactions 

We have searched for evidence of binary interaction during 
X-ray flares by dividing our sample of flares into spectro- 
scopic binaries and all other systems, i.e. wide (or visual) mul- 
tiples and single stars, in which such interactions can not take 
place. The comparison of flare rates F showed that large X- 
ray flares are significantly more frequent on spectroscopic bi- 
naries: FsB = 1.43 ± 0.25% and Fno„-sB = 0.64 ± 0.12%. 
We have taken account of all components in multiple systems 
when evaluating the flare rate. Therefore, the difference in F 
between close binaries and other stars seems indeed to indi- 
cate that magnetic interactions within close binaries leads to 
increased flare activity. But note, that interbinary events are ex- 
pected to have longer durations because of the larger scale of 
the magnetic configuration. Our statistical observations did not 
show an increase of the time scales for spectroscopic binaries. 



9.2.6. Spectral signatures during flares 

From the lower panels of Fig. [ill it can be concluded that for 
most of the observed events the spectral hardness has increased 
during the flare. Due to the large uncertainties, however, the 
changes in the mean hardness are only marginal. But, note, that 
the uncertainties represent the standard deviation (computed by 
taking into account upper/lower limits to the hardness) and thus 
reproduce the spread in the data. 

We think that the X-ray emission of TTSs is harder than that 
of Pleiades and Hyades stars (see Table pi]) for two reasons: (i) 
Because of their circumstellar envelope TTSs suffer from much 
stronger absorption than Pleiads and Hyads, and absorption is 
stronger for Pleiads than for Hyads due to the larger distance of 
the former, (ii) the younger the stars, the stronger the activity, 
and therefore the harder the spectrum. 

10. Conclusions 

We have determined flare rates for PMS stars, Pleiades and 
Hyades on a large data set and found that all stars are observed 
during flares for less than 1% of the observing time. Both fre- 
quency and strength of large X-ray flares decline after the PMS 
phase. 

To probe whether the activity changes in the presence of a 
circumstellar disk, e.g. as a result of magnetic interactions be- 
tween the star and the disk, we have compared flares on cTTSs 
and wTTSs. We find that flares on cTTSs are stronger and more 
frequent. 

A comparison of flares on spectroscopic binaries to flares 
on all other stars of our sample shows that the flare rate is by a 
factor of ^ 2 higher for the close binaries. 

The flare rate of fast rotators is enhanced by a factor of 3 
as compared to slowly rotating stars. 

To summarize, our analysis confirms that age and rotation 
influence the magnetic activity of late-type stars. All previous 
studies in this field have focused on the quiescent X-ray emis- 
sion. Now, for the first time the rotation-activity-age connec- 
tion has been examined for X-ray flares. Furthermore, from 
the sample of flares investigated here we find evidence that 
magnetic activity goes beyond solar-type coronal activity: On 
young stars interactions between the star and a circumstellar 
disk or the magnetic fields of close binary stars may play a 
role. 
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